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We construct an analytical expression for the distribution of gluons in the state immedi- 
ately following a heavy ion collision in the quasi-classical limit of QCD given by McLerran- 
^ ' Venugopalan model. The resulting gluon number distribution function includes the effects 

C^j ■ of all multiple rescatterings of gluons with the nucleons of both colliding nuclei. The typical 

l/") ■ transverse momentum k± of the produced gluons is shown to be of the order of the sat- 

(N ■ uration scale of the nuclei Q s , as predicted by Mueller. We analyze the properties of the 

obtained distribution and demonstrate that due to multiple rescatterings it remains finite 
(up to logarithms of k±) in the soft transverse momentum limit of k± Q s unlike the 
usual perturbative initial conditions given by collinear factorization. We calculate the total 
number of produced gluons and show that it is proportional to the total number of glu- 
ons inside the nuclear wave function before the collision with the proportionality coefficient 
Ph! ess 2 In 2. 

^ . I. INTRODUCTION 

Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory has recently became oper- 
ational with new data already being produced One of the challenges facing theoretical heavy ion 
community is the correct interpretation of the newly obtained data. It has been conjectured that in a high 
energy heavy ion collision a thcrmalized state of quarks and gluons, usually referred to as quark-gluon 
plasma (QGP), may be produced ||, allowing us to explore the properties of QCD matter under extreme 
conditions. Understanding the experimental signatures of QGP is a very difficult task. The first important 
step in that direction is finding a correct description of the state of gluons and quarks immediately following 
a heavy ion collision but preceding the anticipated onset of thermalization. Distribution of partons in this 
state has become known as the initial conditions for QGP formation. The succeeding evolution of this quark- 
gluon system and its possible thermalization can be studied by Monte-Carlo simulations 0), by invoking the 
formalism of transport theory S or by analytical estimates of j(||7]] . 

One may try to model the initial condition by considering pairwise interactions of nucleons in the colliding 
nuclei. The particle production would be described by applying collinear factorization formalism to each 
individual nucleon-nucleon collision and superimposing the results Unfortunately the approach has 
several shortcomings. In the integration over the transverse momentum of one of the produced partons one 
is forced to introduce an infrared cutoff po in order to make the integral finite. The resulting production cross 
section depends very strongly on the cutoff po which makes it difficult to determine the initial conditions 
with high precision. Another (related) problem of the purely collinear factorization approach is that it 
does not include the effects of multiple rescatterings of the produced partons with each other and other 
nucleons, which seem to be, at least intuitively, characteristic to heavy ion collisions where the number of 
scattering particles is large. This phenomenon is intimately related to the problems of nuclear shadowing and 
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higher twist resummation. There have been made several attempts to correct the problem by redefining the 
nucleons' parton distributions to include the nuclear saturation effects and then using collinear factorization 
mechanism to describe particle production |J . 

An interesting model of nuclear collisions has been proposed recently by McLerran and Venugopalan , 
which was designed to both explain the nuclear shadowing phenomenon and to provide us with the strategy 
of deriving the initial state parton distributions free of problems typical to collinear factorization approaches. 
The model is based on the observation that at very high energies the parton densities in large nuclei reach 
saturation and the number of partons becomes very large. Then due to a large number of color charge 
sources the gluon emission can be described by the solution of the classical Yang-Mills equations of motion 




nucleus nucleons 




non-Abelian Weizsacker- Williams 
field 

FIG. 1. Non-Abelian Weizsacker- Williams field of a large nucleus, as derived in |l5|Jl6|Jl^] . 

To describe the saturation philosophy let us start by considering a small- a; gluon distribution function 
of a single ultrarelativistic nucleus. Let us suppose that Bjorken x is small, x <C 1, but not too small for 
quantum corrections to become important ■ The quantum corrections at small x bring in the powers 

of a s In 1 jx and are resummed for the case of one rescattering by the BFKL equation Q] . If a s In 1/x <C 1 
we can describe the gluon distribution in the nuclear wave function including multiple rescatterings on the 
nucleons by the classical gluon field calculated in light cone gauge The field has been found in fl5| . [l(| 
and was dubbed the non-Abelian Weizsacker- Williams field of a large nucleus. In the A + = light cone 
gauge (with the nucleus moving in the "plus" direction along the light cone) the field has only transverse 
components which are given by [nBUlfl] 



A ww (x,X-) = I d 2 z dz-O(z--x-) p a {z,z-) ,~ == Sofe^-lT^Vfez-) (1) 
J \x-z\ 

with the gauge rotation matrix given by a path-ordered integral 

S {x,x-) = Vexp(igT a J d 2 z dz-6(z- - jg_) p a (z, z_) Lnflsc - z\n) \ . (2) 

p a is a color charge density operator normalized according to 

(p a (x, x.) p\y, = ^- p(x, x-) 8(x- - S 2 (x - y) S ab (3) 

with p(x, x-) the normal nuclear density in the infinite momentum frame of the nucleus, obeying 

d 2 x dx- p(x,x-) = A. (4) 



Here we have used the notation of to describe the non-Abelian Weizsacker- Williams field. Feynman 
diagrams corresponding to this classical field were found in [Q and are depicted in Fig. As was discussed 
in some detail in the multiple rescatterings are included in the non-Abelian Weizsacker- Williams field in 
the form of gauge rotations. In Fig. II] the field of the rightmost nucleon effectively gets rotated by the fields of 
all the other nucleons in the nucleus (see Fig. [l]). As was argued in |lq| classical approximation corresponds to 
the limit of no more than two gluons interacting with each nucleon, which formally means resummation of all 
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powers of the parameter a? A 1 / 3 with A the atomic number of the nucleus. Each additional power of a 2 A 1 / 3 
corresponds to an extra rescattering and resummation of all of such terms corresponds to resummation of 
multiple rescatterings. For a large nucleus with a 2 A 1 / 3 ~ 1 all multiple rescatterings are important. 

The unintegrated gluon distribution function of a nucleus is proportional to Fourier transform of the 
correlator of two Weizsacker- Williams fields in the nuclear wave function, which was calculated in |lj , |l7j 
and is given by Eq. ( |2l| ) below. This classical unintegrated gluon distribution has a very peculiar features: 
at large values of transverse momentum it falls off like l/fcj., which is a usual perturbative result. As 
the transverse momentum becomes smaller the distribution increases. However below certain momentum 
scale, which is called the saturation scale and is given below by Eq. (|l9|), the growth of the classical gluon 
distribution slows down to just a logarithmic increase, proportional to In Q s /kj_. This suggests that multiple 
rescatterings may help us to avoid the singularities of collinear factorization by introducing the saturation 
scale Q s , which effectively regulates the parton distributions in the soft momentum region. In terms of this 
saturation scale parameter resummation of multiple rescatterings corresponds to resummation of powers of 
Q^/kj_, since Q 2 ~ a 2 A 1 / 3 (see Eq. (|l9|)). The scale determining the value of the strong coupling constant 
a s is also the saturation momentum Q s . Therefore applicability of perturbative QCD and quasi-classical 
physics depends strongly on how big Q s is. If Q 2 ^> ^qcd tnen a s (Q 2 ) <C 1 and the physics described 
above is applicable to the nuclear scattering process. 

As the energy increases (or, equivalently, as we go towards smaller values of x) the quantum correc- 
tions become important. Multiple rescatterings enhanced by quantum corrections become multiple pomeron 
exchanges. There were developed several techniques which resum multiple pomeron exchanges. The tech- 
niques based on resummation of successive classical emissions in the framework of effective lagrangian of 
McLerran-Venugopalan model led to a renormalization group functional differential equation of |19f| . There 
are different effective lagrangian approaches, such as the one developed by Lipatov and collaborators [2C| 
and another one by Balitsky [ pl| . Finally there is an integral evolution equation which was obtained in [ 23 j 
using the techniques of Mueller's dipole model |22| . The equation is similar to the GLR equation of [11 
and was also obtained in [^lj for the evolution of the Wilson lines' correlators. The approach of also 



seems to be converging to the equation of [[23j (see 24 ). Nevertheless, as was argued in [25|,|l0[ the net effect 



of quantum corrections is to increase the saturation scale Q s without significantly changing the shape and 
main qualitative features of the classical gluon distribution. 




FIG. 2. A collision of two ultrarelativistic nuclei at high energies. 



But what does the saturation physics teach us about the initial conditions for the heavy ion collisions? 
Similarly to how it was done for the case of a single nucleus structure functions one has to start by considering 
purely classical case, neglecting the quantum corrections. This of course corresponds to not extremely high 
energy scattering, a s Ins <C 1. For realistic high energies one has to include quantum evolution corrections, 
but classical initial conditions (in the QCD evolution sense) are necessary in order to do so. Quantum 
corrections may again only modify the saturation scale in the classical distribution without significantly 
changing its shape, as it happened for nuclear structure functions [^5 10 1. The problem of inclusion of 
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quantum corrections still has not been solved. However, since the quantum evolution can be represented as 
a series of classical emissions |l^,|22j the qualitative picture of particle production has been constructed in 
p6f and certain physical consequences, such as multiplicity correlations in rapidity in the produced particle 
spectrum have been predicted [ 26| . 

The classical gluon production problem for nuclear or hadronic collisions has been formulated in pTfl by 
Kovner, McLerran and Weigert. They consider scattering of two ultrarelativistic nuclei from McLerran- 
Venugopalan model (see Fig. ||). The valence quarks of the nuclei just pass through each other during the 
collision without deflection from their straight line light cone trajectories. Corrections to this approximation 
are proportional to positive powers of x, and x -C 1. However, the gluonic field generated by the collision 
has non-zero field strength in the forward light cone and contributes to the gluon production. If we would 
like to obtain an expression for the distribution of produced gluon which includes all powers of a 2 A 1 / 3 (or, 
equivalently, Q s /k±), we can do it by solving classical Yang-Mills equation with the ultrarelativistic nuclei 
providing us with the source current p7[ . The gluon field given by the solution of the classical equations of 
motion in the forward light cone would describe the gluon production (see Fig. ||), and, consequently, the 
initial conditions for heavy ion collisions. This is the problem we are going to address in this paper. 

The classical field of two nuclei in the forward light cone has been found in the usual perturbation theory 
to order g 3 in [p7|-p9|. The answer for the distribution of produced gluons is proportional to QsiQfa/^i.) 
where Q s i and Q S 2 are the saturation scales of the colliding nuclei. This is the first (lowest order) term of 
the expansion in powers of Q 2 sl 2 /k\ and corresponds in this sense to proton-proton scattering. The answer 
agrees with the production rate one would get from employing the so-called Lipatov effective vertex |l4||Tl|l , 
or, equivalently, with the result of Gunion and Bertsch p0| . 

A gluon production cross section for a slightly more complicated case of proton-nucleus interactions was 
derived in (l?]]. In the formal language of the saturation scale parameters Q s i and Q S 2 that cross section 
includes all powers of Q 2 sl jk\ keeping only the leading power of Q 2 2 /k\. 

The problem of finding the classical gluon field in nucleus-nucleus collisions has been addressed in lattice 
simulations of Krasnitz and Venugopalan |3l| ] . The numerical distribution of the produced particles has been 
found and exhibited saturation properties expected, such as finiteness in the small transverse momentum 
limit. Mueller suggested that the total number of the produced gluons has to be proportional to the 
total number of gluons in the nuclear wave function with the proportionality coefficient c of order one. The 
numerical value of the coefficient was determined in j3lj to be c = 1.29 ± 0.09. 

In this paper we are going to write down an analytical expression for the distribution of produced gluons 
which resums all powers of both Q 2 sl jk\ and Q 2 s ijk\. That is we are going to address the problem of 
nucleus- nucleus scattering analytically. The paper is organized as follows: in Sect. II we will formulate 
the problem of classical gluon production and make some useful observations revealing the advantages and 
disadvantages of viewing the process in different gauges. In Sect. IIIA we will review the solution of the 
proton- nucleus (pA) problem in covariant gauge which was given in [[l7), showing how multiple final state 
rescatterings play crucial role in gluon production. In Sect. IIIB we will analyze the same process of gluon 
production in pA collisions in light cone gauge and demonstrate that an entirely different set of initial state 
interactions is important there. We will outline certain important cancellations of diagrams, which would 
allow us to write down an expression for the produced gluons' distribution in nucleus- nucleus (A A) collisions 
in Sect. IV. This distribution is given by Eq. (|35|) and is the central result of this paper. It gives the 
distribution of gluons in the state immediately following a heavy ion collision providing initial condition for 
possible thermalization of the gluonic system at later times. We explore the properties of the distribution 
of Eq. ( |35| ) in Sect. V. There we first obtain a simplified expression for the distribution in the case of not 
very large transverse momenta k±_ < Q s , which is given by Eq. (|38|). We demonstrate that the distribution 
of produced gluons ( |35| ) is finite up to \nQ s /k± in the soft transverse momentum limit and scales as l/kj_ 
when the transverse momentum gets large. We calculate the typical momentum of the produced gluons and 
find that for the case of identical nuclei it is of the order of the saturation scale Q s , as was conjectured by 
Mueller in Finally we estimate the coefficient c from |7| and find c ~ 2 In 2, which is close to the result 
of In Sect. VI we estimate the saturation scale using the new data from PHOBOS experiment at RHIC 
and obtain Q 2 ss 2 GeV 2 , which is marginally in the perturbative QCD region and agrees with the result 
of p2]. We end the paper with a summary of the results obtained. 
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II. FORMULATION OF THE PROBLEM 



In this section we are going to review the formulation of the problem of finding the classical gluon field of 
two colliding nuclei and make some observations which will be useful later. 

As was originally stated in j27]] one needs to solve the classical Yang-Mills equations 

d^ v = r , (5) 

with the current J M arising due to the valence quarks in the colliding nuclei. The valence quarks move ultra- 
relativistically along the straight lines on the light cone and do not get deflected in the collision. (Deflection 
is suppressed by a power of the center of mass energy of the colliding system.) Thus the current generated 
by them has non-zero light cone components J + and J_ and zero transverse component J = 0. In a non- 
Abclian theory the source current is not gauge invariant, it gets rotated under gauge transformations. We 
will first construct the current in a particular gauge — covariant gauge 

9 M = 0. (6) 

Before the collision the nuclei do not see each other and do not interact. The current is given by the sum 
of the currents of two free nuclei on the light cone. As was shown in [ p>|]2£| ] the current of a free nucleus in 
covariant gauge is given by a simple superposition of the currents of the point color charges (valence quarks) 
in the nucleus, each of them being parametrically of order g in strong coupling constant. For example, in 
the model of quarkonium nucleus considered in |jl^,^8|, where the nucleus was envisaged as an ensemble of 
point color charges, with each nucleon consisting of two valence quarks (a quark and an antiquark), the free 
current is 

At 

joov (°) =g J2 T a (T?) [5(x. x,.) 5{x - £,.) - S(x. - x\_) S{x - £)] , (7a) 

A 2 

joov (0) = g ^ T a (fa) [s{x+ _ y . +) §{3 , 6{x+ _ y , + ) 6{3 , _ ^ ; (?b) 

J=l 

J cov (0) = . (7c) 

Here x^s (a^'s) and y^s (y^-'s) are positions of the quarks (antiquarks) in the ith nucleon of the first and 
in the jth nucleon of the second nucleus correspondingly. The first nucleus is moving in the "+" direction, 
while the second is moving in the "— " direction (Fig. |^). (T") and (T") are SU(3) color generators of the 
ith nucleon in the first nucleus and of the jth nucleon in the second nucleus. Ai and A^ are the atomic 
numbers of the nuclei. 

Eq. (^) implies conservation of the classical current J M 

J M = 0. (8) 

GOV 



Similarly to what was done in ]28| we will use this condition to construct the covariant gauge current J' f 
during and after the collision. The valence quarks do not get deflected from their light cone trajectories in 
the collision. The only non-zero components of the current are therefore J+ and </_. We can rewrite Eq. (||) 
as 

d+ J cm + d- J c + ov - ig[A c ° v , J c _ ov ] - ig[A mv , J c + ov ] = 0. (9) 
Eq. (P) can be satisfied with the following ansatz for the current 

Jl ov = U^ix)^ {0 \x)U(x), J mv = S- x {x)JT {0 \x)S{x), (10) 

where 

U(x) = Pexp(-i 5 f + dx' + A c ° v ) (11) 
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and 



S(x) = Pexp -ig / dx'_ AT , (12) 



with A± v the components of the unknown solution of Eq. JET). Thus generally speaking the matrices U(x) 
and S(x) are not known. In arriving at the solution of Eq. (g) given by Eqs. (|l0|), ( pi] ) and (12) we have 
fixed the initial conditions for Eq. (^): we required that before the collision the current J c ° v should be given 

by the free nuclear current J^ ov ^ of Eq. (0). This is just a casuality requirement which makes sure that 
there is no interactions between the nuclei prior to collision. Eqs. (|To|), ( pT| ) and (12) certainly satisfy this 
initial condition, since before the collision U = S = 1, as the fields of free nuclei are non-zero only on 
light cone (see for instance (l^^]). Since Eq. (^) is a linear differential equation the solution of Eq. ( |To| ) is 
unique for the given initial condition of Eq. (|7|) . 

Eq. ( pp| ) has a very simple physical interpretation: the valence quarks in the nuclei during and after the 
collision are still moving along the same straight lines on the light cone. The only effect of the collision on 
these valence quarks is the rotation of their color charges by the gluon field created in the collision and by 
the gluon field of the other nucleus. This has been discussed and illustrated in |^8| at the lowest nontrivial 
order in a s . 

Eqs. (|5|) and (|l^) provide us with complete formulation of the problem in covariant gauge. We have to 
solve the Yang-Mills equations (ra) with the conserved current (|hJ). We are now going to demonstrate an 
interesting property of the current. 

Let us perform a gauge transformation with the matrix S(x). The new gluon field will be given by 

A]f = S ATS" 1 - % - (3 M S) S- 1 . (13) 

As easy to see A+ c — 0, which means that the gauge transformation with the matrix S(x) transforms the 
field into the light cone gauge. The current in the light cone gauge is 

Jl c = S(x) U~\x) J c + ov {0) (x) U(x) S-^x), J L _ C = Jt ov {0 \x). (14) 

From Eq. (|l4|) it follows that the "— " component of the current in the A+ c = light cone gauge remains 
unchanged throughout the collision and is equal to the order g "initial" free nucleus current in the covariant 
gauge. That means that the charges of the second nucleus in the light cone gauge do not get rotated in the 
collision. Let us illustrate what this statement means in terms of diagrams. 




A B C D E 



FIG. 3. Diagrams contributing to the classical gluon field in covariant gauge at order g 3 . 

In [p7|-p9| the classical gluon field of two colliding nuclei was found perturbatively at the lowest non-trivial 
order in g, which happened to be order g 3 . The diagrams contributing to the gluon field at this order in 
covariant gauge are shown in Fig. ^. There we present a collision of two ultrarelativistic nucleons, the upper 
one of which is moving in the light cone "plus" direction while the lower one is moving in the "minus" 
direction. The straight lines in Fig. [?] correspond to the valence quarks. The cross denotes the point in 
coordinate space where one measures the field. 

The interpretation of the diagrams of Fig. || has been given in |28|]. In diagram A two gluon fields merge 
to produce the final field. In the diagrams B and C the current of the upper quark gets rotated by the field 



G 



of the lower quark and a gluon is emitted off the modified current. In the diagrams D and E the opposite 
happens: the current of the lower quark gets rotated by the field of the upper quark and emits a gluon. In 
both cases of B,C and D,E the current of one of the quarks gets a rotational correction of the order g 2 (one 
gluon exchange contribution) which could also be obtained from Eq. ( |l0| ) by perturbative expansion of the 
matrices U(x) and S(x) to the lowest non-trivial order [g8j. 




A B 

FIG. 4. A+ = light cone gauge diagrams: A contributes to the classical gluon field while B does not. 

Now looking at the light cone gauge current given by Eq. (|l4|) we see that J^ c does not get any rotational 
corrections whatsoever and remains at the lowest order in g (free nucleus current). Therefore we can conclude 
that, for instance, diagrams D and E of Fig. || do not contribute to the classical gluon field in A + = light 
cone gauge, which is a self-evident statement. However, we can draw more general and less trivial conclusions 
from Eq. (|14|). Since the current of the second nucleus does not get rotated it means that it only rotates the 
current and/or the field of another nucleus. This also implies that there is no diagrams with more than one 
gluon line connecting to any quark line in the second nucleus (except for virtual diagrams where two gluons 
can connect to a single quark line - see |f8|]). The statement is illustrated in Fig. |[ Diagram in Fig. |]A has 
only one gluon line attaching to the quark line in the second (lower) nucleus moving in the "minus" direction 
and therefore may contribute to the classical field. The diagram of Fig. has two gluons interacting with 
the quark in a nucleon of the second nucleus, and the nucleon remains in a color non-neutral state at the 
end. This diagram does not contribute to the classical field according to what we have shown in Eq. (|l4|) . 



III. PROTON-NUCLEUS COLLISIONS REVISITED 



Before addressing the issue of nucleus-nucleus collisions (AA) let us first consider a somewhat easier 
problem of proton-nucleus collisions (pA). Below we are going to review the solution of the problem in 
covariant gauge given in p"^ ] and then proceed by analyzing the same pA process in the light cone gauge of 
the nucleus. Our interpretation of the underlying light cone gauge physics will be slightly different from the 
one presented in JT^]. 



A. Covariant Gauge 



Consider a collision of an ultrarelativistic nucleus moving in the "plus" light cone direction and a proton 
moving in the "minus" light cone direction. For pA collisions we want to solve the same problem of classical 
gluon production as was stated above for AA collisions. In this subsection we will just follow the discussion 
of fT?!] . We will work in A_ = gauge with polarization vector also taken in that gauge, e_ = 0, which 
for the nucleus moving in the "plus" direction is equivalent to covariant gauge (d ■ A — 0) |l7j] . Following 
fl7f we will consider the process in the rest frame of the nucleus and perform the calculations in the light 
cone perturbation theory (see |54| and references therein). Then the physical picture of the gluon production 
is the following: the incoming proton may already have a gluon in its light cone wave function before the 
collision with the nucleus and the system of the proton and gluon multiply rescatters on the nucleons in 
the nucleus. Alternatively the proton can emit the gluon after the multiple rescatterings in the nucleus. 
The diagrams where the gluon is emitted during the proton's passing through the nucleus are suppressed by 
powers of its large light cone momentum i.e., by powers of center of mass energy of the system (eikonal 
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approximation) |17||. Multiple rescatterings are easier to resum by calculating the amplitude in the transverse 
coordinate space [ |l3|fr^ ] . To obtain the gluon production cross section we have transform the amplitude into 
the momentum space and square it. The diagrams contributing to the gluon production cross section are 
shown in Fig. ||. The graph in Fig. ||A corresponds to the square of the amplitude corresponding to the case 
when the gluon is present in proton's wave function before the collision. The diagram in Fig. gives the 
interference term between the amplitude from Fig. |^A and the amplitude in which the gluon is emitted by 
the proton after the collision. Of course a diagram complex conjugate to Fig. p3 should also be included. 
It can be shown that the square of the diagram with late gluon emission does not have any interactions 
in it and can be neglected. (The gluon exchanges between the proton and the nucleus cancel.) In [|l7| the 
interactions counting was a little different from the one we will present below and that diagram was included, 
leading to the same result. 

A 



nucleons 




proton 



B 




FIG. 5. Covariant gauge (or more accurately A- = gauge) gluon production diagrams for proton-nucleus collision 
as considered in Multiple rescatterings in the nucleus determine the interactions in this gauge. 

Note also that in the quasi-classical approximation depicted in Fig. || the interaction is modeled by single 
and double gluon exchanges. The limit of no more than two gluons per nucleon is imposed Jl8| . If a particular 
nucleon exchanges a gluon with the rest of the system in the amplitude then it has to exchange a gluon in 
the complex conjugate amplitude to remain color neutral. Alternatively the nucleon can exchange two gluon 
in the amplitude (complex conjugate amplitude) , but then it can not interact in the complex conjugate 
amplitude (amplitude). This is done in the spirit of the quasi-classical approximation resumming all powers 
of a 2 s A 1 / 3 , as was discussed in the Introduction. 

In the graph of Fig. HA the nucleons of the nucleus interact with both the proton and the gluon by gluon 
exchanges. It was noticed in Jl7[ | that the interactions with the proton can be neglected due to real-virtual 
cancellation. Moving a gluon exchanged between the nucleus and the proton across the cut does not change 
the momentum of the produced gluon in Fig. but does change the sign of the whole term, causing the 
cancellation. That is why we have to consider only the interactions with the gluon in Fig. |^A. Similar kind of 
cancellation does not happen in Fig. ||B. Moving an exchanged (Coulomb) gluon across the cut would force 
us to move it across the gluon emission vertex for the produced gluon on the right hand side, thus changing 
the momentum of the produced gluon. Thus all the possible interactions have to be included in Fig. pp. On 
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the right hand side of the diagram in Fig. |^B only the interactions with the proton are possible. 

To obtain the answer for the gluon production cross section in pA in the quasi-classical approximation we 
have to convolute the wave function of the proton with a soft gluon in it with the Glauber-type propagator. 
The answer has been derived in [[l7| (see also j}5| and references therein). The diagram in Fig. ||A gives JlTj 

and the diagram in Fig. plus its complex conjugate after a somewhat more sophicticated calculation Jl7j 
gives @ 

= I f d 2 bd 2 xd 2 y 1 e *-<z-vj (l - e ^ 2 ^ + 1 - e~y- ^) . (15b) 

d z k dy n J (2tt) tt xfy z V / 

The total gluon production cross section is equal to the sum of the terms in Eq. ([I 

d** A da\ A da vA 

= 1 =±2-. (16) 

d 2 k dy d 2 k dy d 2 k dy 

In Eq. ( |l5| ) x and y are the transverse coordinates of the gluon in the amplitude and the complex conjugate 
amplitude correspondingly counted with respect to the transverse position of the quark in the proton off 
which the gluon is emitted (0). b is the impact parameter, k is the gluon's transverse momentum. Same as 
in (It]] we use a shorthand notation (see also (k| ) 

x 2 Ql = x 2 8 -^f P xG(x,l/x 2 ), (17) 

with p the density of Eq. (Q) taken in the nuclear rest frame. In the two gluon approximating the gluon 
distribution function of a nucleon is 

xG{x,l/x 2 ) = ^^lnJ-, (18) 

with \i some infrared cutoff. For simplicity of calculations throughout the paper we assume that the nucleus 
has a cylindrical shape, with radius R and the height of the cylinder 2R. The cylinder is lined up along the 
z axis. Generalization of our results to a spherical nucleus is trivial. 

In general the saturation scale Q 2 has to be found from the following implicit equation fl^|l3|| 

2 8TT 2 a s N c R 2 

Qs = n 2 -i pxG ( x >Qs>- ( 19 > 

However, since the logarithm in Eq. ( |l§| ) is a slowly varying function we can assume that in our classical 
approximation without any QCD evolution in the structure functions the gluon distribution function is 
approximately a constant and the right hand side of Eq. ( p"9| ) is independent of Q s . Thus Eq. ( p^ ) turns 
from implicit equation into an equality. 

To summarize the results reviewed here we note that the gluon production in pA collisions in covariant 
gauge is driven by the multiple final state interactions. Now we will explore how this picture changes in the 
light cone gauge of the nucleus. 



B. Light Cone Gauge 

Let us now consider pA scattering in the A + = light cone gauge with polarization vector also taken 
in the same gauge e + = 0. A direct analysis of the light cone gauge diagrams could be a little difficult 



p7[ . We are going to use a different strategy, which was already employed previously in |18|. We know the 
answer for the gluon production cross section given by Eq. (|l5|). Here we are going to guess the diagrams 
in the A + — light cone gauge which give us the same answer. Once we guessed the correct diagrams we 
can conclude that the remaining diagrams should cancel with each other, since they do not contribute to 
the cross section. 
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B 




FIG. 6. Gluon production in proton-nucleus collisions in A+ — light cone gauge (see text). 



The light cone gauge diagrams contributing to the gluon production cross section in proton-nucleus col- 
lisions are depicted in Fig. o. The pA scattering process could be viewed in either rest frame of the proton 
or in the center of mass frame. Again we are going to perform the calculation in the framework of the 
light cone perturbation theory. Similar to the covariant gauge case considered above the incoming nucleus 
can emit a gluon in its wave function either before or after the collision with the proton. The one gluon 
light cone wave function of an ultrarelativistic nucleus is given by A (x) ■ e, with A w the non- Abelian 
Weizsacker- Williams field of the nucleus given by Eq. (0) with suppressed X- dependence fl5|,|l6| and e the 
polarization vector in light cone gauge. (One can see that it is the case for instance by calculating the one 
gluon rescattering diagram on the left hand side of Fig. |^A with this wave function and by using conventional 
perturbation theory. Both results give the same answer.) Diagrammatically the light cone wave function 
corresponds to the same set of diagrams as was depicted above in Fig. [|. The fields of the nucleons in the 
nucleus "gauge rotate" the Weizsiicker- Williams field of one of the nucleons jl8| . The interaction with the 
proton can only be by the means of single or double gluon exchanges, as was shown in Sect. II. Eq. ( |l4| ) 
shows that the "minus" component of the current does not get rotated implying that there could not be 
more than two gluons exchanged with the proton in the A + = light cone gauge. 

Before colliding with the proton the nucleus can develop the Weizsacker- Williams one gluon light cone 
wave function which then interacts with the proton by means of one or two gluon exchanges, according to the 
rules of the quasi-classical approximation p8| , |l7f . The square of the graph corresponding to this scenario is 
shown in Fig. ^A. As in Fig. ||the interactions of the proton with the nucleons in the nucleus cancel through 
the real-virtual cancellation leaving only the interactions with the gluon line. One may notice that the final 
state interactions are left out in the diagram of Fig. |A, but as we will show below we do not need them to 
reproduce the contribution of the graph in Fig. ||A, which implies that they cancel with each other. 

The second possible scenario corresponds to the case when there is no gluon in the nuclear wave function 
by the time the collision happens and the gluon is emitted by the nucleus after the interaction with the 
proton. Then the nuclear wave function without an emitted gluon corresponds to the fields of the nucleons 
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rotating the current of one of the nucleons in the nucleus. This is shown on the right hand side of Fig. ^B. 
The nucleon then interacts with the proton by exchanging one or two gluons with it. After that the nucleus 
can emit a gluon to be produced in the final state. Another possibility which is not shown in Fig. but 
which contributes to the gluon production corresponds to the case when the Weizsacker- Williams gluon is 
present in the nuclear wave function by the time of the collision, similar to Fig. [j|A, but after the interaction 
with the proton the gluon merges into the quark line of one of the nucleons, which later re-emits the gluon. 
We could not find an a priori argument prohibiting an emission of the whole Weizsacker- Williams field after 
the interaction. However, as we will see below one needs to emit only one gluon to be able to reproduce the 
results of the previous section. The square of the diagram on the right hand side of Fig. ||B is zero since 
the interactions cancel due to real-virtual cancellation |L7|]. The only contribution we get from it is the 
interference term depicted in Fig. ^B. There on the left hand side we have the same diagram as in Fig. |^A 
except that now interactions of the proton with the "last" nucleon in the nucleus do not cancel. We will 
show that the diagram of Fig. ^3 provides us with the contribution equal to that of the graph in Fig. ||B. 

Let us now calculate the diagrams in Fig. ||. The contribution of Fig. ||A can be obtained by convoluting 
the correlation function of the fields on both sides of the cut with the gluon-proton interactions amplitude. 
The result yields 

da LC\ _ f & x d2 V ik-(x-y) 2 / ,w w . AWWf, A \ -u s n 2 N c _ 2 ^ 



, , /(/ ~ j r2 , a ' e - J Tr (A ww (x) • A ww {y)) fe - yf xG{x t 1/fe - yf). (20) 



In [|16|J17| 1 the correlation function of two non-Abelian Weizsacker- Williams fields in the nuclear wave function 
was found to be 

Tr(A ww (x)-A ww (y)) = ^— - (l - e -fe-*W*y (21) 

N — ' TTa s (x — y) z \ J 

Employing Eqs. ( ^l|) and (|l^) in Eq. (|2(i| ) and defining new variables z — x — y and b = y we obtain 

a& ■ Z**^^^ (***•-»)• < 22 » 

Using Eqs. (63) and (64) from jlTj one can see that 

i l f 2 v ■ (z + y) 



In -^-z = - Id y =y — — =j (23) 
zfix 1 it J yf{z + yy 

where the y integration is cut off by 1/ fx in the infrared limit. Inserting Eq. ( |23| ) into Eq. ( p2| ) and comparing 
the result to Eq. dlqa) one can see that 



dcr LCi _ d<7 i A 



d?k dy d 2 k dy 



(24) 



Thus we have shown that the contribution of the diagrams in Fig. ^|A is equal to the contribution of the 
diagrams in Fig. [l5|A. 

The calculation of the graphs depicted in Fig. is a little more complicated. Similar to diagrams of 
Fig. |^A a correlator of two gluonic fields is involved. However, in the field on the right hand side of the cut 
in Fig. ||b only rotations of the source are allowed. Similarly to what was done in jl8|,[l7j we argue that the 
fields of the nucleons rotate the quark line to which the emitted final state gluon is attached, as well as the 
Coulomb gluon's field coming from the proton below. Everything is gauge rotated, except for the gluon line 
of the emitted gluon. After applying Ward identities we conclude that effectively only the vertex where the 
emitted gluon connects to the quark line on the right hand side of Fig. [|B is rotated. This corresponds to 
rotation of the source of the gluon field emitted [H . The effect of this gauge rotation is modification of the 
expression for the field from Eq. |lj) into 

A££d(x,x-) = [ d 2 z dz-0(z--x-) p a (z,z^) .-~-. 7 So(z,z^)T a S Q 1 (z,z^). (25) 
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After we sum over the all possible connections of the Coulomb gluon lines to the "last" nucleon in the 
nucleus in Fig. |^B the resulting expression would depend on the transverse coordinate of the quark in that 
nucleon to which the Coulomb gluon couples. Thus we can not factorize the averaging in the nuclear wave 
function from the interaction terms anymore, like it was done in obtaining Eq. (GO). With all the above- 
mentioned complications in mind we write the contribution of the diagram in Fig. |gBas (see for a similar 
calculation) 

d 2 kdy tt J (2tt) 2 J - \\x~z\ 2 \y-z'\ 2p U ' JP { ~' -> 

xTr [S (x, z.)T a S^{x, Z-)S ^, z'_)T h (z' , z'_)] [x - z) 2 xG(x, l/(x - z) 2 )"j (26) 

Similar to what was done before in Jl5|,[l7j we argue that we can average the densities p a (z_, Z-) p b (z_' , z'_) in 
the nuclear wave function independently and we can employ Eq. (|^) to do so. In |17| the following result 
was derived (see Eq. (48) in p7| ) 

(Tr [S a {x,z-)T a S^{x,z^)S a {z,z^T a S^{z,z^)\) = 
= C F N C exp (-^y (x - z) 2 p(x, zS) xG{x, l/(x - z) 2 ) (z_ + , (27) 



with ±zo- is the upper (lower) limit of the z_ integration in Eq. (|26|). In Eq. (27) we assume that p(x, zJ) 
does not vary much between x and z_, which is justified for a large nucleus. Using Eqs. (||) and (|27j) in 
Eq. (Eq) we obtain 



d <C2 = I f d' 2 x d 2 y e ik.(x-y) d 2 z <XsC F x-z y-z , _ Q , /4 > ^ 

d 2 k dy 7r J (27r) 2 7r \x— z\ 2 \y — z| 2 V 



Defining new variables x = x — z, y — y — z and b — z and adding the contribution of the complex conjugate 
to Fig. ||B diagram (da 1 ^ 3 /d 2 k dy) we can compare the result with Eq. (|l5|b) and conclude that 

7 pA j pA rlrr pA 

aa LC2 _j_ aa LC3 _ " g 2+3 ^29) 

d 2 k dy d 2 k dy d 2 k dy 

We have thus shown that the contribution of the diagrams in Fig. ||B is equal to the contribution of the 
diagrams in Fig. ||B. 

We have proved that the diagrams in Fig. ^ are the only diagrams contributing to the gluon production 
cross section in A + = light cone gauge. Thus a whole class of diagrams with final state interactions does 
not contribute to the cross section. Several examples of such graphs are shown in Fig. ^. 

The diagram in Fig. represents a class of diagrams where the non-Abelian Weizsacker- Williams wave 
function of the nucleus after interaction with the proton merges with another (non-interacting) Weizsacker- 
Williams wave function. The graph in Fig. can be viewed as a similar to Fig. (?]A process, where after 
proton-nucleon interaction a gluon field is emitted, which later on merges with the non-Abelian Weizsacker- 
Williams wave function of the nucleus. There is another class of diagrams which do not contribute to the 
cross section where two Weizsacker- Williams wave functions merge with each other and with a Coulomb 
gluon coming from the proton through a four-gluon vertex, producing a gluon in the final state. Those 
diagrams are probably suppressed in the old-fashioned light cone perturbation theory as requiring the gluon 
fields' merger to happen at a particular light cone time when the system is passing the proton. 

As was shown above all of the final state interactions shown in Fig. ^ do not contribute to the gluon 
production cross section in pA collisions. Therefore they should cancel, either with each other or individually 
due to some other cancellation mechanism. From considering gluon production in pA collisions in the light 
cone gauge we may draw the following conclusion: the gluons produced by the collision do not merge with 
the non-Abelian Weizsacker- Williams wave function of the nucleus. 
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FIG. 7. Some of the final state interaction diagrams which do not contribute to gluon production in pA collisions 
in A+ = light cone gauge. 

Another important diagram which a priori should be contributing to the gluon production in pA scattering 
is shown in Fig. |[ The diagram contains virtual interaction between the proton and the field of the nucleus. 
Even though there is no interaction on the right hand side of the diagram it is allowed by the ruled of the 
old fashioned light cone perturbation theory, where the energy is not conserved in the vertices [fl7|,j34|| . The 
diagram of Fig. || has a remarkable feature in it, which was absent in the graphs shown in Fig. fn merger 
of two produced gluons. In Fig. ^ two Weizsacker- Williams wave functions of the nucleus first interact with 
the proton producing gluons, which then merge with each other. This merging is different from the ones 
considered in Fig. ^. There one of the merging gluons was produced in the interaction with the proton while 
the other one was just given by the non-interacting Weizsacker- Williams wave function of the nucleus. In 
Fig. H both merging gluons were first "produced" by the interactions and then merged together. Here again, 
since the diagram of Fig. ^| does not contribute to the gluon production cross section it has to either be zero 
or cancel with some other diagrams. This gives us a very strong reason to conclude that the gluons produced 
during the collision do not merge with each other at later times in A+ = light cone gauge. Even a more 
general conclusion can be conjectured: gluons produced in the interactions do not interact with any other 
gluons afterwards in A + = gauge. 



FIG. 8. A virtual diagram which does not contribute to the gluon production process in pA collisions in A+ — 
light cone gauge. Triple gluon vertices are marked with black dots. 

To conclude we note that the physical interpretation of the scattering process appears to be gauge depen- 
dent: in the covariant gauge case considered in the previous section the gluon production was dominated 
by multiple final state rescatterings. In the light cone gauge multiple rescatterings vanish. The information 
about them is now contained in the light cone wave function of the nucleus. The same observation about 
the interplay of initial and final state interactions was made in Jl7[ for the case of current-nucleus scattering 
with the current j = — \F^. 



Now we can employ the results we have obtained in Sect. Ill to write down an ansatz for the distribution 
of produced gluons in nucleus-nucleus collisions. Let us consider a head-on central collision of two ultrarel- 
ativistic nuclei, as was shown in Fig. 0. We will be working in A + = light cone gauge. The calculations 



nucleus 




proton 



IV. NUCLEUS-NUCLEUS COLLISIONS 
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can be done either in the center of mass frame or in the rest frame of one of the nuclei. We will work in 
the rest frame of the second nucleus. The diagrams contributing to gluon production in AA are depicted 
in Fig. ||. They are somewhat similar to the diagrams of Fig. |^ and of Fig. |5[ The incoming nucleus may 
or may not have a Weizsacker- Williams gluon in it. In the first case the system multiply rescatters in the 
second nucleus at rest. This is illustrated in Fig. |]A. Similar to pA case multiple rescatterings between the 
first and the second nuclei cancel. Only the interactions with the gluon survive. The interference graph of 
the amplitude from Fig. ^A and the amplitude where the gluon is emitted after the interaction is shown 
in Fig. ^B. Analogous to pA case we work with diagrams where the final state interactions are limited to 
multiple rescatterings in the second nucleus and a single gluon emission (or absorption) by the first nucleus. 
As was demonstrated in Sect. Ill in the case of proton- nucleus scattering all other final state interactions 
including "produced" gluons' merging cancel. Here we argue that this also happens in the nucleus- nucleus 
collisions. 

A 

nucleus #1 




nucleus #2 



nucleus #1 




FIG. 9. Diagrams contributing to the gluon production in nucleus-nucleus collisions in the A + — light cone 
gauge. 

In Sect. Ill, analyzing proton-nucleus scattering we concluded that the diagrams where the gluon produced 
by interaction with the proton merges with the non- Abelian Weizsacker- Williams wave function of the nucleus 
cancel. That allowed us to neglect similar diagrams in Fig. || above. However, there exists also a somewhat 
different class of diagrams, one of which is shown in Fig. |l0|. There the gluons that merge in the later stages 
of the collision both were produced in the interaction of the Weizsacker- Williams wave function with the 
second nucleus. This brings us back to the class of diagrams depicted in Fig. |^, where there is also a merger 
of two "produced" gluons. Since the diagrams of Fig. ^ canceled in the case of pA collisions we have a strong 
reason to believe that the graphs of the type shown in Fig. |l^ also cancel in AA collisions. Thus Fig. || 
contains all the diagrams that contribute to gluon production in nucleus- nucleus collisions. 



=Q 



=<2 



FIG. 10. An example of a diagram which is not included in the gluon production mechanism of Fig. []. To avoid 
possible confusion the three gluon vertices in this diagram are marked by black dots. 
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We can calculate the diagrams in Fig. ^ along the lines outlined in the previous section. The contribution 
of the graphs of Fig. ^]A is obtained by a simple eikonalization of Eq. (p(f), and, correspondingly Fig. |]A. 
The result for the number distribution of the produced gluons yields 



dN l A = f d2x d 2 V ik-(x-y) 2 rjY / jWW foA \ f e -(x-y) 2 Q 2 s2 /-l 

cPkdy J (27r) n 



- Tr {AY W (x) ■ AY W (y)) 1 U-^-?) Q ^ 4 - 1 . (30) 



Here the indices 1 and 2 denote the fields and the saturation scales of the first and the second nuclei 
correspondingly as well as over which nucleus wave function the quantities are being averaged. With the 
help of Eq. (EHI) we rewrite Eq. (5o|) as 



dN i = 2C F f d xdy ^(x-y) 1 A _ e -(x-yf QUA ( x _ e -(x-y) 2 Q 2 s2 /A . (31 ) 
d 2 kdy a s TT 2 J (27r) 2 \x — y\ 2 V /V / 

As was demonstrated in Appendix A of |l7) multiple rescatterings of Fig. also just exponentiate. There- 
fore, in order to calculate the contribution of the diagram of Fig. ||B we should substitute 

|^fe-z) 2 xG(x,l/(x-z) 2 ) (32) 

in Eq. © by 

(33) 



Performing all the same integrations and averaging that were done in obtaining Eq. (28) we end up with 

dN f A — l^E. f d ' x ^ Jkte-v) d 2 b^.lL 1 ( i _ e -£ Q 2 jA (i _ e -2 2 (M ) 
d 2 kdy ~ a s TT 2 J (2tt)3 ~ x 2 f^^^V ) V )' {6 *> 

To get the full answer we have to add the contribution of the diagram which is complex conjugate to the 
one shown in Fig. |^B. Summing up Eqs. (|3l|), ( |34| ) and its complex conjugate, we obtain the expression for 
the number distribution of gluons produced in a heavy ion collision 

dN AA = f d?z ik . z 1 / -J qua 

d?kd?bdy a s Tr 2 \ J (2tt) 2 z 2 



d 2 x d 2 y ik .(z:-y) jL_ 
(2tt) 3 6 x 2 



y~ X" m r^- 



V 



In 



l 



-V 2 



)( 



1 - e 



-v 2 Q 2 2 /4 



(35) 



Eq. (35) is our main result. It provides us with the number of gluons produced in a head-on zero impact 
parameter heavy ion collision per unit transverse momentum phase space, per unit rapidity interval at the 
given impact parameter b. Eq. ( |35"| ) includes the nucleon-nucleon scattering result of Fig. ||| calculated in 
|p7| f29f as well as the proton-nucleus scattering contribution of Fig. ^ found in JItJ . We are going to explore 
the properties of the distribution (B5J) in the next section. 



V. PROPERTIES OF THE CLASSICAL DISTRIBUTION 

Let us evaluate Eq. ([55]) in the approximation in which we neglect all the logarithms of the transverse 
coordinates |l7[] , since logarithm is a slowly varying function and can be assumed to be a constant compared 
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to powers. That technically means putting In -A- 



lnr4- 



1 in Eq. (|35|). That also concerns terms like 

x 2 Q 2 , which in general also have logarithms of |sc| in them, as follows from Eqs. (1^) and There we 

also put the logarithms to be of the order of one, similar to how it was done in p7 1 . We have to note that 
this approximation is good only for not very large transverse momenta k±<Q s . When gluon's momentum is 
large, k± 3> Q s , the logarithms of the transverse coordinate are crucial for deriving the correct asymptotics 
of the distribution function of Eq. (|3F 
Employing the fact that 



d 2 y 
(2tt) 2 



-ik-y y_ 



i k 



(36) 



in the second term of Eq. (|35|), integrating over the angles of z and x and performing similar integrations in 
the third term of Eq. j35|) we get 



dN 



AA 



c F r 00 dx 

d 2 kd 2 bdy a s ir 3 J Q x 



J: 



Integrating over x in Eq. (37) we find 
dN AA C F 1 



d 2 kd 2 bdy a s 27r 3 k 2 



_ k 

12 „~ 



Qtie 



(37) 



(38) 



The distribution of Eq. (|38|) is plotted in Fig. |llj as a function of k/Q s for the case of two identical cylindrical 
nuclei with Q s \ = Q S 2 = Q s and with the cross sectional area S± = irR 2 w 50 fm 2 w 1250 GeV -2 . Note 
again that Eq. (|3q) is valid only in the not very large transverse momentum region k±<Q s . 
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FIG. 11. Distribution of the produced gluons given by Eq. ( |38[ ) for a central A A collision as a function of k/Q a 
with the transverse area of the nuclei being S± = 50 fm 2 . The approximation of Eq. feq) is valid only for k± Q s . 



As one can see the distribution in Fig. |Tl| remains finite as k±/Q s 
of Eq. (pq) then the distribution goes to a constant 



0. If one takes the k± <C Q s limit 



dN 



AA 



C, 



d 2 kd 2 bdy 



a s 27r 3 



as 



oT <<L 



(39) 



This means that the exact expression of Eq. ( p5| ) may only have logarithmic divergences in the infrared 
limit. This conclusion is very interesting, since the initial conditions for heavy ion collisions which are 
generated by pairwise interactions between nucleons in nuclei without multiple rescatterings have a power- 
law divergences in the infrared limit [g . Even the pA gluon production cross section of Eq. ([l5|) diverges as 
& t small transverse momenta because it includes multiple rescatterings in only one nucleus, since one 
nucleus is involved in the scattering process. Therefore finiteness of Eq. (B8h at small transverse momenta 
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demonstrates that multiple rescatterings are the reason the hadronic and nuclear single particle inclusive 
production cross sections remain finite in the soft momentum region. 

Since the distribution of Eq. contains the lowest order in a s diagrams in it (see Fig. |J) one readily 
derives that in the k±/Q s — > oo limit the distribution falls off as l/k\ p7|-p9t 



dN AA Q 2 sl Q 2 2 k_ 

d?kd?bdy ~ a s k 4 ^ Q 



as -f- -> oo (40) 



which is a well-known perturbative result. As one can see extrapolation of the usual perturbative expression 
of Eq. ( |40|) into the soft momentum region would lead to singularities and strong cutoff dependence of the 
total number of the produced gluons. The multiple rescatterings of Eq. ( |35| ) resolve this problem. 

A simple calculation shows that the typical transverse momentum of the gluons in the distribution of 
Eq. (|38f) is given by 

(k 2 ) = 2 QliQ 2 s2 (41) 

{ ~ ' (Qii + <&) HQli + <&) - Q 2 si mQ^ - Q 2 2 lnQ 2 s2 ' 
For two identical cylindrical nuclei Eq. ( pilf ) gives 

<i*> - & (42) 

That is, the typical transverse momentum of the produced gluons is of the order of the saturation scale, as 
was conjectured by Mueller in (7). Since the saturation scale for a large nucleus scales as Q 2 ~ A 1 / 3, with 



atomic number, as could be seen for instance from Eq. (19), it may get quite large, much larger than the 
non-perturbative QCD scale Aqcd- Then most of the produced gluons would have momenta high above 
Aqcd which would justify the use of perturbative QCD in the problem [^0|J^]. 

Finally, the total number of gluons produced in the collision can be found by integrating Eq. ( |38|) over 
k±. The result yields 

= ^2 [(<& + QI2) HQli + QI2) - Qli - Q 2 2 ^2] , (43) 

which, for the case of identical nuclei gives 

dN AA C F Ql In 2 
d 2 bdy a s TT 2 



(44) 



In Mueller suggested that in a high energy nuclear collision the gluons in the wave function of the incident 
nucleus get liberated by the interactions with the nucleus at rest. The coherence of the incoming nucleus 
gluonic wave function is broken by the second nucleus. Thus the total number of produced gluons should 
be proportional to the total number of gluons in the wave function of one of the nuclei before the collision 
with the proportionality coefficient c, which should be of order one 0. Therefore one may write 

which, using Eq. (|2l]) can be rewritten as 

^ c (46) 



d 2 bdy a s 2iT 



Comparing Eq. (46) to Eq. (|44j) we conclude that 

c = 2 In 2 w 1.39, (47) 

which is very close to the result of the numerical estimates of Krasnitz and Venugopalan giving c = 1.29 ± 
0.09 |3^]. The obtained value for the "gluon liberation" coefficient c = 2 In 2 is close to one, as was originally 
suggested by Mueller Q . 
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VI. DISCUSSION 



Eq. ( ]44| ) allows us to estimate the saturation scale Q s knowing the multiplicity of produced particles. Of 
course one should be careful in interpreting this estimate, since at high energies the purely classical picture 
considered here breaks down and quantum corrections bringing in powers of a s lns become important. 
It has been conjectured though (h]^(| that these corrections would not change Eq. (Q) and would only 
(considerably) increase the value of the saturation scale Q s on the right hand side of it. Another issue one 
should be worried about in this kind of an estimate is that the classical gluon production picture presented 
here does not include the interactions at late times, which may lead to thermalization of quarks and gluons 
produced, and may also modify the total number of gluons |6| . Eq. gives us the total number of gluons 
immediately after the collision, which may be different from what the detectors count at the end due to the 
importance of 2 — > 3 and 3 — > 2 processes at the later stages of the collision || . Also Eq. (Q) gives us the 
total number of gluons produced, which is not quite equal to the total number of pions, kaons and other 
hadrons observed in the detector. Here we will just assume that due to entropy conservation the numbers are 
very close to each other p6[ . Keeping all the above mentioned restrictions in mind we may nevertheless try 
to estimate Q s in our classical picture here using the newly emerging RHIC data [Q. PHOBOS experiment 
has measured total charge multiplicity per unit pseudorapidity in Au+Au collisions yielding the result [jjj 

j Au+Au 

— ^ = 555 ± 12(stat) ± 35(syst) (48) 

at the center of mass energy y/s = 130 AGeV. In our crude estimate we will multiply the number given 
in Eq. (|48|) by 3/2 to account for charge neutral particles and use the resulting number as a lower bound 
estimate of dN/dy in Eq. (E3). (dN/dy is a little larger than dN/dr].) Again we use a cylindrical nucleus 
approximation with the cross sectional area S± = 50 fm 2 . We assume that the strong coupling constant is 
a s ~ 0.3. (In general a s = a s (Q 2 s ) and we have to treat Eq. (Q) as an implicit equation.) The result for 
the saturation scale is 

Ql w 2.1 GcV 2 for Au + Au at s/s = 130 AGeV, (49) 

which is close to and even a little larger than the estimate of [§2] . The saturation scale of Eq. (E^) appears 
to be marginally in the perturbative region. As energy of the RHIC beam reaches 200 AGeV the particle 
multiplicity will increase too, leading to an even larger saturation scale, which would make the use of 
perturbative QCD at RHIC even better justified. 

To summarize the results of this paper we repeat again that we have derived the classical distribution 
of gluons produced in the ultrarelativistic heavy ion collision (Eq. (|35|)), thus constructing classical initial 
conditions for the evolution of the gluon system leading to a possible gluon thermalization. It would be 
very interesting and important to analyze the subsequent evolution of the gluonic system in the framework 
of McLerran-Venugopalan model and see whether the onset of thermalization is possible before the system 
falls apart and to what experimental consequences that would lead. Important first steps in that direction 
have already been made [|5|-fjJ. Eq. (|35l) can also be applied to describe minijet production in the proton- 
proton collisions at very high energies | p7p9| ]. The proton's high energy wave function consists of many sea 
partons which may serve as color charge sources for the classical field similar to nucleons in the nuclear case 
p0| , ^9|j2l| . We have derived a simplified expression for the distribution of Eq. ( |35| ) which is given in Eq . ( |38| ) . 
We have demonstrated that the distribution is finite in the soft transverse momentum region (Eq . (|39|)) and 
approaches the usual perturbative result when the transverse momentum becomes large (Eq. (Eo|)) . We have 
shown that the typical transverse momentum of the produced gluons in the collision of two identical nuclei 
is of the order of the saturation scale Q 8 (Eq. ([42])). Finally we have observed that the total number of the 
produced gluons is proportional to the total number of the gluons in the nuclear wave function with the 
proportionality coefficient c = 2 In 2 (Eq. (|47|)). 
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